On the suppression of the diffusion and the 
quantum nature of a cavity mode. Optical 
bistability; forces and friction in driven cavities 

Karim Murr§ 

INFM, Dipartimento di Fisica, Universita di Camerino, 1-62032 Camerino, Italy 

Abstract. A new analytical method is presented here, offering a physical view 
of driven cavities where the external field cannot be neglected. We introduce a 
new dimensionless complex parameter, intrinsically linked to the cooperativity 
parameter of optical bistability, and analogous to the scaled Rabbi frequency 
for driven systems where the field is classical. Classes of steady states are 
iteratively constructed and expressions for the diffusion and friction coefficients 
at lowest order also derived. They have in most cases the same mathematical 
form as their free-space analog. The method offers a semiclassical explanation 
for two recent experiments of one atom trapping in a high Q cavity where the 
excited state is significantly saturated. Our results refute both claims of atom 
trapping by a quantized cavity mode, single or not. Finally, it is argued that 
the parameter newly constructed, as well as the groundwork of this method, are 
at least companions of the cooperativity parameter and its mother theory. In 
particular, we lay the stress on the apparently more fundamental role of our 
structure parameter. 
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1. Introduction 

In 1994, Kimble (1994) quoted the necessity of theoretical work to be done in 
order to understand the dynamics of an atom in a cavity where the mode is not 
a prescribed quantity but rather a fully playing actor. A forefront difficulty is to 
draw an efficient theoretical framework that treats the strong coupling regime in 
open systems, externally fed, and dissipating by cavity decay k and/or spontaneous 
emission 7. The first area of cavity QED is usually devoted to coherent aspects of 
the coupling g between the atom and the cavity mode. It is well treated within the 
dressed states formalism (Cohen- Tannoudji 1990), a vision that contains in block 
all the information about the quantized field plus the atom and where, at first, 
dissipation and driving are not accounted for. On the other hand, driving a cavity 
leads to dissipative structures. Those are well understood within the semiclassical 
theory of optical bistability (Lugiato and Narducci 1990). Semiclassical means here 
that the field is factorized from the internal state of the atom, in the operatorial 
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sense. This corresponds in bistability theory to treat the field classically by means 
of Maxwell-Bloch equations. The smallness of the decay rates (k, 7) allows the field 
to be kept inside the resonator boundaries, and hence, a long interactivity occurs 
before the escape of one photon or one atom. A difficulty arises as there is then 
a competition between the scaled Rabbi frequency of the dressed states (coherent 
part) and the Bonifacio-Lugiato cooperativity parameter C = g 2 /2nj (Bonifacio and 
Lugiato 1978,1975). 

Recently, two major experiments reported the trapping of an atom in a cavity 
containing about one photon on average (Hood et al 2000, Pinkse et al 2000). In 
both papers, the trapping times of about a millisecond and the registered signals 
were considered to be a signature of the quantum nature of the field. Another 
study (Doherty et al 2000) suggested by theoretical and numerical arguments that the 
experiment of (Pinkse et al 2000) is well understood without quantization assumption 
of the field. One argument was to compare the diffusion coefficient and the trapping 
potential which they found to be close to the free space ones (Gordon and Ashkin 
1980). In contrast, the experiment of (Hood et al 2000) deals with a parameter 
regime where the diffusion is suppressed by nearly a factor of ten when compared to 
the free space case. (Hood et al 2000, Doherty et al 2000) conclude that such high 
suppression is a signature of the dynamics embedded in the Jaynes-Cummings ladder 
of dressed states, pointing out indirectly on the quantum nature of the trapping field. 
(Hood et al 2000) support their conclusions by citing other experiments, including 
(Hood et al 1998), that found good agreement between the full quantum calculus and 
the experimental data of a heterodyne signal, on the single atom level, whereas the 
semiclassical bistability state equation fails. 

In parallel, (Pinkse et al 2000) attribute the trapping mechanisms to cavity- 
induced cooling (Horak et al 1997, Hechenblaikner et al 1998). That mechanism is 
based on the weak-field limit and is extendable (Domokos et al 2001) to higher photon 
number provided a low saturation is kept. Their work follows and precedes several 
studies concerned with the mechanical effects of light on atoms within a cavity (see 
e.g Doherty et al 1997, Vuletic and Chu 2000). A better understanding together with 
an extension of the analytical calculus to higher intensities and/or saturation, on the 
single atom level, seems difficult as one expects, either the inclusion of multiphoton 
processes within the dressed states formalism, or an extension of the bistability state 
equation. 

Having drawn the different areas to be studied, we bear in mind that the need 
for the full quantum calculation to fit some experimental data does not necessarily 
mean that one actually deals with quantum mechanics, rather, it indicates that the 
known semiclassical theories are insufficient. We demonstrate in this paper that both 
experiments of (Pinkse et al 2000, Hood et al 2000) are substantially explained by a 
semiclassical approach. A strong result is that, even if one tunes the external laser 
resonantly to a dressed state, neither experiment needs the ladder of dressed states 
to be explained. That is true for elementary steady state quantities as well as for 
the diffusion and friction coefficients. To that purpose, we develop a fully analytical 
framework that allows one to construct the needed quantities, structure and dynamics, 
by taking increasing orders of the coupling <?, without any dressed states consideration. 
We introduce a new physical parameter, until now un-noticed, which is efficient for 
the treatment of driven dissipative systems, and by essence is a cooperative quantity. 
The method, in this paper, is limited to the explicit construction of factorized states 
that are an iterative version of the bistability state equation. That is possible by 
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introducing the idea of the referred states. The referred states allow the construction 
of the steady states. Six steady states are derived, their common structure is given 
in section [3 Those are then grouped into two distinct families, the bounced states 
section |31 and the polarized states section 0] typically valid for increasing saturation 
and photon number. Our parameter is introduced in subsection 11.21 and is discussed 
in subsection l2.3l Section [3 draws mathematical remarks on the role of our parameter 
in bistability theory. Finally, we derive three expressions for the diffusion in sectional 
and two expressions for the friction in sectional Diffusion and friction are derived in 
such a way as to be independent of any explicit expression of the states. 

1.1. The physical system 

In dealing with cavity QED one usually uses the driven Jaynes-Cummings Hamiltonian 
for a 2 states atom represented by the lowering and rising operators (<7,cp), a =\g)(e\, 
where \g) and |e) stand for the ground and excited states respectively, and a single 
mode with annihilation and creation operators (a, a^). We consider the so-called 
rotating wave approximation and focus only on the internal state of the atom 
interacting with one cavity mode. The master equation for the reduced density matrix 
p of the system is governed by the Liouvillian C (h = 1): 

p = Cp, Cp= -i[H,p]+"/L a p + nL a p , (la) 

with the dissipative parts written in the form L p = 2bpw — pb^b — b^bp, for any 
operator b, and with the Hamiltonian (in the frame rotating with the probe frequency 

^probe) • 

H = oj a cr j o- + ujca^a + g{a)a + acr t ) + Ea) + E*a . (16) 

Here oj a and lu c are the atomic and cavity frequencies respectively, detuned with 
respect to the probe frequency u) pro b e - Although these detunings are usually labelled 
by (A a ,A c ), we keep our notations in order to have the formulas more compact: 
(u> a = uj atom— u probe, u c = ui cav — Lo pr obe) ■ 7 and k represent respectively the atomic 
and cavity-mode decay rates, E is the strength of the probe field that coherently 
drives the cavity mode, and g is the coupling that we assume real and dependent on 
the position of the center of mass of the atom. For illustration purposes we assume a 
cosine dependence along the cavity axis x, g(x) — gocos(2Trx/X) 1 with wavelength A and 
go is the maximum of g. All the results presented are independent of any particular 
form of the coupling g. Finally, the dipole force operator Tdip reads: 

Fdip = -{^9)Fd, F d = aJ+cta. (lc) 

1.2. Relevant dimensionless physical parameters 

In order to reach a consistent understanding of the physics of the system, we look for 
the relevant parameters to be formed from (|la|lljlfejl . For E = 0, a natural description is 
provided by the ladder of dressed states, whose structure is basically understood by the 
square of the scaled Rabbi frequency = g 2 (n+1) / S 2 = C n ( Cohen- Tannoudji 1990), 
where S = uj a — lu c is the detuning mismatch and n the mode quantum number. C n is 
independent of the probe frequency w pro be and built by atomic and mode parameters, 
where the latter has a discrete quantum nature. If one includes a driving, the steady 
state mean photon number in the cavity without the atom is given by the complex 
amplitude «o : 

A^o=|ao| 2 , a = E/uj c , (uj c = uj c - in) . (2) 
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The parameter ao depends upon the probe through the strength E and the cavity 
detuning lo c . When the driving is sufficiently high E < g , a given dressed state 
becomes coupled to upper and lower manifolds. It becomes hard to even deduce the 
mean photon number Nq from this formalism. Indeed, ao is a continuous description 
of the cavity mode that tends to dismantle (dilute) the discrete dressed states by 
connecting them. When the atom is driven by a classical field —a (e.g — ao), with a 
strength — ga, the relevant quantity is g/uj a (u> a — u a — ij) where \ui a \ 2 /g 2 gives the 
amount of (normalized) field intensity needed to saturate the atom. Symmetrically, 
when the atom is described by a classical dipole oscillator with polarization a — > /3, of 
strength g/3, a relevant parameter for the cavity mode is g/ui c , whose square modulus 
represents the amount of field intensity that is exchanged between the cavity and the 
atom (compare with ao ©)• 

The parameters {E /ui Cl g /ui a , g /ui c ) depend on the probe frequency w pro t>e and, 
contrary to C„, they do not describe any settled structure of the coupled-atom-mode 
system. A parameter accounting for that is the cooperativity parameter C = g 2 /2k^ 
in bistability theory, but this quantity does not include the detunings. So we use the 
ratios g/to a and g/ui c to introduce a new dimensionless complex parameter: 



The parameter v presents the same structure as the dressed states parameter C n except 
that it depends on the probe frequency oj pro be and contains no quantum signature of 
the cavity mode. We shall stress throughout this paper that, when the cavity mode 
is described classically, v has a similar status as the (squared) scaled Rabbi frequency 
C n . As v depends on the probe frequency ujprobe, then it is likely to be related to 
the driving strength E. We show this below. Several properties for v should be 
mentioned: First it can be seen as the ratio (g/u c )/(u} a /g), when squared it compares 
the number of photons (average) provided by the atom to the number of photons 
needed to saturate the atom. In dynamical terms, it measures the ratio between 
the atomic dispersive shift g 2 /u> a and the field complex frequency uj c . v reduces to 
the cooperativity parameter —2C when the atom and cavity are at resonance with 
the probe, and finally its imaginary part contains the mode-pulling formula in laser 
theories Q(^) oc nuj a + jlu c (Lugiato and Narducci 1990). Finally, the deep meaning 
for v is mathematically dealt with in subsection 12.31 it constitutes in this paper a 
central physical parameter. 

1.3. A problem with the standard semiclassical factorization procedure 

A description of the system that avoids the dressed states formalism could be provided 
by the bistability state equation, whose underlying theory is generally understood by 
invoking the so-called semiclassical approximation. This approximation assumes a 
factorization of the mean values of the product of atomic operators A and mode 
operators C. It is satisfied as far as the state of the system can be factorized 
into a product p = p m ® p a , where (p m , p a ) are states for the mode (m) and the 
atom (a) respectively. With the notations, tr(CAp) =< CA >, tr(Cp m ) =< C >, 
tr(Ap a ) =<A>, that reads: 



<CA>=<C><A> 



(4) 
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The Heisenberg equations of motion for (<a> 1 <d->,<a z >), with the population 
difference a z = a^a — aa^ >=<&>* is implicit), are written: 

<&>= — iuj a (<a> — — <aa z >) (5a) 
<a>~ — iuj c (<a> +a + — <a>) (56) 

<& z >= -27 - 27 <a z > +2ig(<a)cr> - <acr t >) . (5c) 

By the use of the closed system is solved at steady state; by eliminating the 
atomic variables (< a >, < >, < a z >), one obtains the bistability state equation for 
the output field a Q b = — <a>: 

v g 2 
aob = ct Q {l-— s ob — 2— —\a b\ 2 , (6) 

l + s ob \UJ a \ Z 

where s b stands for the corresponding atomic saturation parameter, <<r z >= —1/(1+ 
s b). Notice that, written in terms of v (|2J), the equation obtained is much clearer 
in form than the usual expression (Lugiato and Narducci 1990, also see section Q. 
That equation is cubic in the steady state output intensity |< a >\ 2 = |a f,| 2 , thus 
it would possibly lead to a hysteresis cycle. The non-linear S-shape of the bistability 
state equation is up to now unobserved on the single atom level. (Hood et al 1998) 
demonstrated a clear disagreement between the measured heterodyne detection and 
the semiclassical bistability state equation, while the data are still well explained by 
the exact steady state of the system ifTall fTTit . 

2. General formalism 

2.1. Iterative factorization procedure: Definition of the referred states 

Our criticism to the above procedure is based on the remark that the expectation 
values (< a>,<cr>,<cr^ >, < a z >) are all considered as variables in (|5a|) (|5^j) (|5cf) . 
Whether they appear directly or they are issued from the factorization assumption 

Mainly there is no other choice than to deduce the explicit values by solving the 
system to obtain equation JJjJ. 

Our approach is still to construct factorized states, but we treat each dynamical 

equation I|5a|llj56|l l |5c|l independently and iteratively, with assumption (0J viewed 

differently for each operator product of the form CA, at each step. A specific 

example is provided to illustrate our method. A low saturation regime is assumed, 

the exact steady state p s corresponds to an atom almost in the ground state \g), 

< a z > s « —1 (vanishing saturation parameter). Moreover, the exact mean values 

of product operators CA can be written in the form of a product of mean values. 

Eventually, the low saturation regime is chosen such that one has for the bistability 

state equation s b 3> 1, hence © is a bad (factorized state) approximation. Given 

these conditions, we approximate the exact steady state p s by p s ~ p m <8> |<?}(g|, for the 

mean value < aa z > only at this step, < aa z >« — < a >, to obtain < a > (I5ajl . that 

is then replaced in l|56|l to obtain <a>= — ao/(l — v). We then assume that <a)o> 

factorizes according to to deduce <cr z > from (|Fc|) . Those first steps give, 

ao ga . . 

a= , <(7 z >=—l + s, <er>= — , (7a) 

l—V LU a 

where we have defined a = — < a >, and s = 2g 2 /\u> a \ 2 x \a\ 2 is the saturation 
parameter. The value of a just deduced corresponds to the mean value of the field 
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in the very low saturation regime (Hechenblaikner et al 1998). An iteration is then 
started. These equations are valid for s <C 1, so one can put < a z >« —1/(1 + s). 
we then deduce < a > by setting one more time < a >= l\'xil and use Q where 
the product is now treated as < aa z >—< a >< a z >= a/(l + s). Those procedures 
correctly normalize the populations, and we check that (jFc|) is exactly cancelled if we 
use the values (<cr>, <a z >) just deduced and assume (@J for <aV>. The iteration 
is stopped here without treating (|5fcjl a second time, otherwise, as shown in the next 
section, another state is obtained. The final equations are: 

a o 1 get 1 

<a z >=-—— , <°>=^T—- ( 7b ) 



1 — v " 1 + s u> a 1 + s 

Those equations cannot be deduced by (|HJ, where one would set s b = 0, because we 
are assuming s i, 3> 1 (an example is illustrated in the next section). 

At this stage, three remarks are given. Firstly, as we treated equation 1)5 only 
once, the value <a>— —a is the same in (I7a|ll)7&)) and it has been deduced by using 
the ground state \g). As we assume low saturation, \g) is close to the final state 
verifying 1)7 fcj) . but we distinguish it from the latter because strictly speaking one has 
<u z >=/= —1. The ground state \g) plays the role of a referred state that was used to 
construct the final state. The status of the referred states is strengthened below by 
the construction of other states. The second point is that, the atomic mean values 
in 1)76)1 correspond to optical Bloch equations where the atom is driven by a classical 
field of amplitude —a, hence the properties are well known. Finally, the value of a 
determines entirely the atomic quantities, thus all those steps can be reduced to two 
steps: Firstly, one treats (|5a)ll5fe)l by invoking |@J and uses the referred state to deduce 
a — — <a>. From now on we write only the necessary equations: 

— / 9 s 

<a>= —ioj a \<o-> — — <aa z >) 

W ° g ( g ) 
< d >= — ito c {< a > + Q'o + — < o >) . 

Secondly, a Bloch equation description is assumed to deduce the atomic mean values. 
In the next section, an optical Bloch Liouvillian with a damped cavity mode is derived, 
its steady state p a is a factorized state, it is known provided that one gives a. 



2.2. General form for the steady states 

The method is based on a shift of the cavity mode by a position-dependent complex 
number a = a(x): 

a = c — a . (9a) 

Such general shift is then used to describe equivalently the Liouvillian (|la|l with the 
Hamiltonian, 

H = LJ a a* a — gacr^ — got* a + g(c) a+ca^) +uj c c^ c + (E—uJcCt)^ + (E — uj c a)*c , (gty 

with the dissipative terms (kL c , jL a ). Four components of the cavity mode dissipative 
part L a where already included in the Hamiltonian through the imaginary part of the 
complex frequency uj c — u) c — in. If a given a is found such that the contribution to the 
steady state p s from the coupling g(c^ a+ca^) cancels approximately the contribution 



Driven cavities in the strong coupling regime: Structure and dynamics 



7 



from the field driving term (E — uj c a)c ! + (E — ui c a)*c, then p s can be approximated 
by the factorized steady state p a of the Liouvillian C a , 

C-aPu = 0, £ a p a = -i[H a , p a ] + nL c p a + jL a p a = 

t t + (^c) 

H a = uj c c'c + u> a o-'a — gota' — get o . 

The system ()9c|) describes an atom driven by a classical field of amplitude —a, i.e 
optical Bloch equations, and a cavity mode representing a damped displaced field 
(c = — ius c c). The steady state of the cavity mode part of l|9c[l . as much as the 
corresponding diffusion and friction coefficients ||, can be described classically by 
the correspondance operator— ^complex number, c — > —a c (t), a c = —iuj c a c . The 
operatorial description of the field is nevertheless maintained for convenience's sake 
because in (|9d) the atom is quantized. The steady state of C a is a factorized state 
and written p a = Q a |0) (0| where |0) is the vacuum of the displaced cavity mode, 
c |0) = 0, it is a position dependent coherent state of the annihilation operator a with 
the eigenvalue -a, a |0) = —a(x) |0). Q a represents the steady state of the atom and 
is entirely defined by (3{x), which is in turn determined by a(x). One can show that 
(see e.g Cohen- Tannoudji 1990): 

Q a = [W + & + /3* a + \/3\ 2 ] /(l + s), (9d) 

= —a, s = 2|/3| 2 , (C n = w a - i-y) , (9e) 

where s = s(x) is the atomic saturation parameter and l/(l + s) ensures a unit trace. 
The expectation value of an operator O in the steady state p a is given by the trace 
tv(Op a ) =<0> a - Some expectation values to be used in this paper are: 

< a^a > a — \o<-(x)\ 2 , < & z > a z 



1 + s(x) 



(9i) 



where the force was calculated by (|f c(l . Since a state p Q is entirely determined by a 
complex number a = — <a> a we shall refer to a as a state. 

The approximation p s w p Q is valid under the semiclassical approximation. It is in 
general satisfied if the time scale that characterizes the internal evolution of the atom 
is well separated from that of the cavity mode (bad or good cavity limits, large atomic 
or cavity detuning). In practice, we estimate the validity of our states by writing the 
dynamical equations (C a p a = 0): 

< c > a = —iuJ c (a — a + — <a> a ) , <<r> Q =0, <<r z >a=0 (10a) 
<aa z > a = -iuj c {<o z > a (a -a)-— <o> a )) , (106) 

and look under which parameter regime the equation for <c> Q (and <ab~ z > a ) are 
totally or partially cancelled. However, as higher order equations need to be analysed, 
we compare our states p a with the numerical results. 

The determination of the different values for a is grounded on a manipulation 
of the Hamiltonian i|96t (and the dissipative parts). Other values are deduced by 
applying a similar reasoning to the Heisenberg equations of motion JSJ). 

| For those dynamical quantities a similar algebraical Liouvillian is used in sections 



Driven cavities in the strong coupling regime: Structure and dynamics 



8 



2.3. An atom- cavity-mode cooperative picture: (v,E) and the dressed states 

In this subsection, we manipulate the Liouvillian C and show how v appears 
mathematically. We then show a relationship between v and E; and, finally, we relate 
v to the dressed states. The whole reasoning comes with an intuitive picture, as well 
as indications on how states will be built in the next sections. 

v can appear by two consecutive operators shifts. The first shift corresponds to 
the choice a — E/uj c = ao in (|9fc|) . The driving is entirely transferred to the atom, 
— gaoa^ — ga^a, we cancel this term by a second shift, on the atomic operators, 

a = S + , (11) 

with the convenient choice (3(x) — gao/uj a = (3$ (see H9e[) ). The atomic driving is now 
exactly cancelled, and a new term from the coupling Hamiltonian comes out. Such 
mode-driving term, gfiS + g/3*c, is then written in terms of v (j3Jl and E, to give: 

H = cj Q E f S+ ui c Jc + g(JE + cE f ) + vEc* + v*E*c , (kL c , 7 L S ) . (12) 

In order to include the first reaction of the atom on the cavity mode, a third shift 
is operated, c — > c — veto, to cancel the mode driving term vE in (|12fl (result not 
written). Once that done, one returns from £ to a, and checks that these three local 
shifts are equivalent to one global shift, starting from the initial Hamiltonian (|9b|) . 
with a = ao(l + v). As we shall see, this value of a gives one simple and practical 
state. Another cycle will generate terms in v 2 in l|12f> . and so on. By those three local 
shifts, the strengh E has been 'transported' and 'crossed' twice through the coupling 
term in l|9&|) (under braced term). Actually, the coupling term in l|16(l plays the role 
of a pivot, it can be seen as the net in a Ping-Pong game between the atom and the 
cavity mode, where E is the ball and v is the dimensionless quantity that transports 
the ball. 

An implicit relationship between v and a finite value of E is addressed by the 
following two statements. Firstly, for E ^ 0, v has a natural role, as shown in the 
scheme above or in the bistability state equation @). Secondly, this is not the case 
when E = 0, unless for the trivial value v = 1. The second statement is explained 
by setting E — in (|9^j) and by shifting the atomic operators using (fill) . One 
then equates each of the new driving terms, (u5 a /3 — ga)a^ and (g(3 — uj c o)c^ (and 
hermitian conjugates), to be zero^f. By that, it is easy to show that the condition 
a 7^ implies oj a oJ c = g 2 . But, this equality is impossible to satisfy unless one has 
simultaneously k = 7 = (9(f) = 0) and io a Lo c = g 2 (3?(f) = 1). From that brief 
analysis, we extract two physical consequences. Firstly, by E = and K = 7 = 0, 
one has the dressed Hamiltonian ill /tt , and, the condition u> a u> c = g 2 means that one 
is tuning the lowest dressed state |— ) into resonance, in such case one has in fact 
cj_ - ujprobe = (w a + w c)/2 - [4g 2 + (uj a - lo c ) 2 ] 1 ^ 2 /2 — 0. Secondly, in O, as long 
as one decay rate only (k or 7) has a finite value, then v^l and the atom (or mode) 
can never return the same probe strength E to the mode (or atom) (vE 7^ E). As 
E and v are somehow linked, we say that v is also the ball, it is the ball that moves, 
and any exponent of v says how many exchanges took place. We call our method the 
Ping-Pong. Graphically one can see the local shifts as a ping-pong exchange: 

e — » y Ct ° v \, /* ao " 2 \ /* aol/!> (13a) 

A) Pov 0O" 2 

% A second shift is necessary to cancel the Active driving created by the first shift. 
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or the global shifts with the reactions on the field only and where the atom is 
represented by the hooked arrow, 

a a (l + v) <-^> ao(l + u + v 2 ) ... , (136) 

the same last operation symmetric for the atom is of course possible by substituting 
ao for /3q. The starting in 1)13 ojl is explicitly shown on the field, this would be reversed 
if the atom was driven first. 



3. The bounced states 

3.1. First bounced state 

If one performs an infinite iteration in (|13b|) . a known serie in v (complex number) 
leads to the first bounced state a — an,: 

axb = • (14a) 

1 — v 

Those infinite iterations could be avoided by making directly a shift of a in 1)9 fejl 
and by cancelling the atom and mode driving terms. The two conditions obtained 
(E — uj c a + 3/3 =0, u) a (3 — ga) give back (|14al) if uj c w a j^g 2 . 

That state is our example (I7^|l in subsection 12.11 where the derivation has 
been opposed to the one that leads to the bistability state equation ©. We add 
some important comparison with the derivation proposed by (Hechenblaikner et al 
1998). Their calculations purport that, under weak driving, the Hilbert space can be 
restricted to the three states 0, bare) , \g, 1, bare) , |e, 0, bare)}, where 'bare' means 
the Fock states in the basis of a^a. Thus, in the equation for < & > = (JSJ), the 
equality <aa z >= — <a> is strictly satisfied, and all expectation values factorize. By 
<a>=0 ©, and by the replacement <a>— >— au, we obtain l|14a|) . Actually, their 
atomic steady state quantities differ from ours i|14afl l|9c|) owing to the presence of the 
saturation parameter sib = s^0 in the denominators of 

A part from that, firstly and by the operations above, the polarization < a > in 
the equation for < a >= (|HJ is proportional to the field amplitude, < a >oc au, 
and hence it has been replaced by the variable au- Consequently, by looking at the 
equation for < a >= (JHJl, the cavity mode 'sees' itself because it is 'bounced' by 
the atom. We call (|14afl a bounced state. Secondly, such state is known to be valid 
for higher photon number, provided that low saturation is ensured (Domokos et al 
2001). That means that a low saturation regime can be chosen, but with high photon 
number, Ao > 1 and 1 — v\ 2 < 1, in order to have the state almost orthogonal to the 
subspace {\g, 0, bare) , \g, 1, bare) , e, 0, bare)}. Their derivation is restricted to weak 
driving, and hence we question back the operation < aa z >= — < a >, that is an 
exact one in this subspace. The problem is absent in our description, because, firstly, 
the state a±b has its field part described by a coherent state, which is, in the bare 
basis, centered around the Fock state corresponding to a quantum number n« |aih| 2 ; 
secondly, the use of the referred state |<?) gives < aa z >= — < a > without any weak 
driving assumption. 

In figure^ we plot the heterodyne transmission and the probability of the excited 
state as a function of the atomic detuning tv a . The parameters are chosen in order 

+ The authors introduce the saturation parameter, but only briefly in the bad cavity limit, and finally 
linearize their equations by putting s = in the denominators in order to obtain the friction and the 
diffusion coefficients. 
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Figure 1. Comparing the bistability state equation with the exact result and our 
approximations. Variation with respect to the atomic detuning u) a . Figure(a), 
heterodyne transmission. From bottom to top at ui a = go- state !14nl . exact state 
(solid line), state IT^cl . bistability state equationlcTt. Figure(b), probability to be 
in the excited state. Same hierarchy. The parameters are (7 = 0.02,k = 0.6,uj c = 
0.1, 30 = 1,^0 = 0.37). 



to make the bistability state equation no longer valid (s b 3> 1), whereas our state 
(and the one derived below) reproduce reasonably well the exact curve. For those 
parameters, one has < aa z > s ^< a > ob < o z > b and < a > t< cr z >ob^ z — < a ><,(,, 
whereas the exact steady state verifies <aa z > s ~ — <a>ifc. For larger detunings, all 
curves meet, and the bistability state equation converges faster than (|14afl . 

Finally, we check the validity of that state. By I jlOaj l. < c >q= e with e = 
iuj c aoi'[l—i']~ 1 x sif,/(l + sib). Thus, e ~ for (su — 0) « or a <C 1 or \v\ <C 1, 
and, by generalization to other Heisenberg equations, p s ~p a . Notice that l|10&(l is 
totally cancelled for the first bounced state, this is also true for <acrt > a ^0. Finally, 
bear in mind the important condition \v\ 3>1, or |1— v\ 2 3>1 (for the photon number), 
that should ensure low photon number and low saturation regimes; aib — > (drop in 
transmission), thus by (|i^ejl [3\b^0 =^-si6— >0. 



3.2. Second bounced state 

This state is constructed by using the first one. As far as low saturation is achieved, 
<|14all works for probabilities of the excited state around <CT t cr> 16 < 0.1. The increase 
of the driving increases the difference Sit>— 0, which we interprete as a departure of 
that state (atomic part) from its own referred state, i.e the ground state \g). The idea 
is, if sit>— is not too large, then one can iterate the procedure by viewing the first 
bounced state as a referred one. Therefore, we write in JSJ < aa z >~< a >< a z >u, 
with the population difference (|14a|) (|9^): 

1 . , 2g 2 (x) |a | 2 

< o-z >i 6 = -— , 8 lb (x) = ,iy . ' ' |2 , 146 

1 + s u \uj a \ 2 \1 - v[x)\ I 

and, as before, one deduces from (JSJ) the second state <a>^ —a.2b, 

a 2 b = , where vi b {x) = — . (14c) 

1-Vlb l + Sib(X) 
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This state is still a bounced state, and we can see that it can be deduced by a ping- 
pong scheme (|13 fo|> with the simple replacement v — * v^. The population difference 
of the referred state an, is now transported. Contrary to the first bounced state, the 
second one tends to «o for increasing driving. The parameter regime where such state 
is valid can be understood by estimating as before (|10all . By (|14c|)(|iTe| )(|9fll, 

S 2b =<cr>2b= — = — — , (14d) 

1 + S 2 b w a l + S 2 b 

with obvious notations, < with s ocao(s2b — su)- Basically such state is valid 

as far as the population of the excited state is close to that of l|14all . In practice, 
it means that by keeping the atomic detuning sufficiently large, the driving can be 
increased to reach < er^er >2b^5 0.2. The intermediate photon number regime can be 
reached. In figure ^ this state is presented and shows better agreement than (|14all 
because the exact population of the excited state corresponds here to values around 

<CT t CT> s «0.15. 



4. The polarized states 

An increase of the population of the excited state leads to the appearance of the two 
level structure of the atom, as seen in the equation < cr^a > s — < cr > s < a > s = 
m(s)/2, where m(s) = s 2 /(l + s) 2 . For low saturation, m(s) ~ s 2 , and hence the 
two level structure of the atom can be represented by the saturation parameter only. 
Bounced states have the property to be defined only by the saturation parameter of 
the referred state, and the increase of the population of the excited state transfers the 
steady state p s from one bounced state to another. In such regimes, the cavity mode 
amplitude is entirely returned by the quasi-pointlike atom (< aa z >—< a>< a z > p ), 
i.e < a > is a variable and p is a fixed (referred) state. Else, the polarization of an 
atom behaves as < a >^ s 1 / 2 for low saturation and as < a >^ s~ x / 2 for very high 
one. Therefore, the trick is to say that a clear appearance of the two level structure 
should be interpreted by treating the polarization <a> in the equation for <a>= 
(jBJ) as one object that the field 'sees' as a whole. In this case the polarization is the 
referred quantity. 



4-1- First polarized state 

The first polarized state corresponds to the first hooked arrow in (|13fr|l : 

ai p = a (l + v) . (15a) 

The reactive term va,® can be written g /uj c x S re f, where S re f is a referred polarization 
of the atom because it contains no information on the actual value of the mode, but 
rather, it is incremented to the initial one cto (S re f = g/u) a cto). The regimes where 
Ijl5all is valid is roughly estimated by calculating <er>i p , to find < c >o= e where 
eoc— h>ao[l— (l+^)/(l+Si p )]. Thus e « for v w s lp . Such a condition is a restriction 
but it can apply to high saturation regimes. 

For the experiments of (Hood et al 2000) and (Pinkse et al 2000), the probe is 
tuned to the lowest dressed state — ) that is, as we said, the particular condition v f« 1. 
As also swl, Ijl5a|l should work to give <a^ a> s «4A^o at an antinode. In figure[21 we 
plot the mean photon number in the steady state as a function of the axial position 
of the atom. The agreement is satisfactory, in particular the maximum and minimum 
that give good agreement with the presented signals in both papers. For (Hood et al 
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Figure 2. Mean photon number in the steady state as a function of the axial 
position x. Figure(a), from (Pinkse et al 2000) (7 = 0.187,/c = 0.087,^ = 2.8,u) c = 
0.31, 90 = 1,^0 = 0.9). Figure(b), from (Hood et al 2000) (7 = 0.02, k = 0.13,<J a = 
1. 13, uj c = 0.7, go = l,JVo = 0.32). The doted lines are |oip| 2 , the analytical formula 
115 al . Solid lines are the exact numerical calculations. 

2000), the transmission reaches about 1.2 — 1.3, while for (Pinkse et al 2000) it goes 
around 4. State l|15afl is not an approximation of l|14a|l . both states meet for v—>0. 

4-2. Second and third polarized states 

More general states can be deduced by iterating the procedures above. Bounced states 
are sensitive to resonances, such as v — *■ 1, whereas the first polarized state (|15 of) is 
not valid for large values of v. The idea is to use both advantages. 

The second polarized state is calculated by reference to the polarization (|14d|l of 
the second bounced state (|14cfl . Consequently, this state is referred in a hierarchy to 
the three states (114 cjl (|14afl and the ground state. With the use of l|14djl . and then 
H14cjl . this gives the second polarized state (<a>= and <a>^ — a 2p ): 

a 2p = a + -^-S 2b = Q!26(l + v q ) v q = vib - v\b , (15&) 

where v 2 b = + S2b)- The last equation in ljl56|) shows that such state is 

also deduced by a ping-pong exchange, with an input field being now a 2 b and a v 
that transports a difference of population difference (compare with (|15afl(|13bjl ). For 
s 2 b — > sib, the state reduces to a 2 b- In general, (|15bjl is more flexible on the detunings, 
it would work for probabilities of the excited state up to 0.25 — 0.3. For situations 
where the first bounced state l|14afl fails, the resulting polarized state l|156|) would 
smooth the results. That is the case for the experiments we discuss. 

When such polarized state does not cancel resonant points, we can use the first 
polarized state (|15af) as the original referred state. To that end, l|15all and |9/|| are used 
to form the saturation parameter Si p of the first polarized state; it is then regarded as 
a referred state to deduce a third bounced state a^b- We then calculate its polarization 
"5*36 |9i|, and view it as a referred state. One obtains thus the third polarized state 
(and the third bounced state): 

a 3p =a + -^-S 3b = "36(1 + v q ) v q = v 3h - v\ P (15c) 




Figure 3. Figure(a), heterodyne signal as a function of the driving strength. 
From (Hood et al 1998), (7 = 0.02, n = 0.33,w a = u> c = 0.166,9 = go = 1). 
From left to right, at ordinate 0.1: bistability state equationJIQ, exact state 
(solid line), state Il5fel . The lowest horizontal line is the first bounced state I14al . 
Figure(b), mean photon number as a function of the axial position. Parameters, 
( 7 = 0.187 , K = 0.094, u) a = 6,oj c = 1/6, g = 1 ,N = 2.63) . From bottom to top: 
statc !15dl . exact state (solid line), state J15 61 . state llBal , Notice the difference 
between the maximum and minimum for the mean photon number (intermediate 
regime). 



a 3 b = —^ — , vi p = — - — , (sip = ^-fi\aa\ 2 \l + u\ 2 ) , (15rf) 

l-Vlp l + Slp \u)aV 

where v%b = v/{l+s^b). The third bounced state l |15t^ . that comes with the derivation, 
could work for the low saturation regimes where the other bounced states do not 
(typically v k, 1). The third polarized state works well for the experiments discussed 
here, in particular for the heterodyne signal. However, as the atomic detuning is low, 
some difficulties are encountered with the regime of (Hood et al 2000) . Some quantum 
effects of the cavity mode are in fact present but we demonstrate below that they are 
not needed for the dynamical quantities. Eventually, those third states are like 115of) . 
as they are not valid for large values of v. 

In figure Ufa) , we plot the heterodyne transmission as a function of the driving 
strength and compare it with three approximations, among which (|15b(l and the 
bistability state equation ©. This plot is figure 5 of (Hood et al 1998), where 
the experimental data fit in well with the theoretical quantum model (|lb|) and it 
constituted a part of the experimental results that supported the conclusion on the 
quantum nature of the field. But equation j|156fl has a classical-field interpretation. 
We mention, however, that other steady state quantities are not well represented for 
these parameters. In figureEHb), we illustrate the intermediate photon number regime, 
and with 3i(i/) < 1. Such parameter regime could be interesting experimentally, as it 
should offer a high signal to noise ratio. The reason why these polarized states work 
efficiently comes from the fact that we applied two different elementary procedures, 
the bounce and the polarization. Therefore, such states are almost stable when the 
parameters are varied. In counterpart, poorer results than with the other states are 
to be expected though. 

In the next two sections are provided several expressions for the diffusion and 
the friction coefficients. It is believed that the arguments presented are easy to 
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understand provided that one pays particular attention to where stationary quantities 
are assumed, in particular the different en's and /3's to be used. We repeat that the 
results derived below are independent of any explicit expression of the steady state 
quantities. 

5. Diffusion at zero velocity 

In the papers of (Hood et al 2000, Doherty et al 2000), it is suggested that the high 
suppression of the diffusion, when compared to the semiclassical free space case, is a 
relevant signature of the quantum nature of the field. As the mean photon number 
is around one, it was concluded that one observed single atoms bound in orbit by 
single photons (Hood et al 2000). The semiclassical free space diffusion coefficient is 
calculated by assuming an equivalent standing wave of strength < a > s , driving the 
atom by g<a> s . The corresponding Liouvillian (Gordon and Ashkin 1980, Cohen- 
Tannoudji 1990) is the atomic part in (|!Tc|) . with —a — ><<z> s . It leads to a diffusion 
about ten times (factor~ 7.5) bigger than the one obtained from the full calculus 
l(Ta |l fT&|l . Consequently, the validity of (|9"cjl is questioned back, for the calculation 
of dynamical quantities, within the cavity QED setting. One answer is addressed 
by pointing at a distinction between dimensionless quantities and dynamical ones. 
Suppose a given Liouvillian C with steady state p 8 , Cp s = 0, with eigenvectors and 
eigenvalues (p k ,E k ), Cp k — E k p k . If it is possible to derive a scaled Liouvillian, 
C\ = XC, then both Liouvillians have the same steady state and eigenvectors, but the 
eigenvalues differ E k = XE k . A dimensionless quantity, such as < a > s , is therefore 
unchanged, while a dimensioned one, as the diffusion, will change according to its 
relation with the scaling procedure. More in accordance with our case is to suppose 
a Liouvillian C = L\ + £2 where (£i,C2) have the same steady state p s . If each 
sub-Liouvillian is scaled by different factors, then C\ = X1C1 +X2C2 still has the same 
steady state p s as (£, £1, £2) but the eigenvectors will this time differ from those of C. 
Dimensionless quantities within the Hilbert space also change, except those that are 
evaluated at the steady state. Starting from the full Liouvillian 1161) . the approximated 
one (|^ej) offers a a good steady state structure, but fails as far as the dynamics are 
concerned. 

5.1. Diffusion by eliminating the mode variables 

A solution is found by first noting that the atomic part of <j!Jc|) can be written by 
using Hll|) . (u) a Y^Y> ,7-ks) for the Hamiltonian and the dissipative part respectively. 
The atomic steady state Q a verifies independently [S^S, Q a ] — 0, and L^Q a = 0. 
An effective atomic Liouvillian is then derived by treating JSJ) in the operatorial form 
(without noise terms) , by setting a = in JSJ and by eliminating a in the equation for 
a. As the equation for & z is unconsidered, then o~ z should be handled carefully. As 
long as o~ z is alone, we keep it as an operator, and when it is in product with a, it is 
'frozen' by assuming that it takes a steady state value —1/(1 + s), given for example 
by some a and |9j. That reads: 

(T = -iil a o- - iga a z , £l a = w a (l - —7—) , 
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which could be derived by the following Liouvillian, 

ffj>*> = n a (aO£tE, r a (x)L s , h a = n a ~ir a 

g 2 (x) 1 



n (x) = w a -t 

r o (ar) = 7 + k- 



C u 2 + K 2 l + s{x) (16a) 
g 2 (x) 1 



K 2 1 + s(x) 

where we already shifted a by an amount flu = gcto/^a, that is not related to 
a by (|9e|l but directly by a supplementary infinite ping-pong exchange l|136jl with 
v — > f/(l + s). As a rule, such a difference will not affect the result (because (3h ~ /?), 
this supplementary bounce procedure can be seen as a way to seek fluctuations around 
the steady state of (|16af) . which is close to the one defined by a. However, if the free 
space steady state of (|9c|l is needed, the procedure will then be stopped one step 
before, and hence a z will be frozen at some value "before the steady state". Actually, 
this procedure is performed so as to avoid non-linear equations. What is obtained is a 
new Liouvillian that is scaled with respect to the free space case, where the dynamical 
parameters (u a ,"f) of the atom are modified by the inclusion of the reaction of the 
field. As we saw, the dynamics is likely to change, but with the same steady state. 
Finally, the dipole force is defined by: 

Fdi P = -(Vff)Frf, F d = -a d a f - a*fj , (166) 

with a dimensionless atomic force operator Fd approximated by assuming a classical 
field characterized by ad , to be precised below. The calculation is easily deduced 
from the free space case. The diffusion tensor Dij is written: 

D * = (^(Vff) x D , 

where D (hereafter called diffusion) is given in |Appendix A.l| The result is close to the 
free space expression, except in a global factor slightly modified, plus a supplementary 
part that is in general negligible. The free-space-like term D at (ctd, Sh) <|A.3|1 can, in 
general, be approximated by D at (a,s) = D at (a) . It indicates that the force l|16fcj) is 
defined by = a, and the saturation parameter is displaced "one step behind", thus 
returned to the steady state value. In figure^a) we plot D at (a, s) for the experiment 
of (Hood et al 2000). The free space case reaches a value around 12 at the antinode 
(not shown). The suppression is well explained. By keeping the term the 
diffusion gets closer to the exact one and even closer if Sh is left as it originally was 
(the diffusion will reach the value 2 at the antinode). 

For comparison's sake, it is also tempting to express the force (|16fr(l in terms of 
ao, since everything is already included in the modified atomic frequency and decay 
rate, and the atom is driven by —gag. In such a case, D at (ao,s) is lower than the 
exact value. Actually, as (|16afl was used, we did not find any proper argument to 
legitimate one value of a d rather than another in the expression of the force l|166j) . 
That is seen in figure ^a) where different values of a d — (a, ao) approach the exact 
result differently. For these experiments, the dominating term in D at (ad, Sh) is the 
known free space diverging term in s^: r^lay^s^l + s^] _1 - Here, formally, it does 
not diverge but reaches a constant value r" 1 ^^ 2 (unless for a d — a). The origin 
of that divergence has been interpreted (Gordon and Ashkin 1980) in terms of the 
dressed state approach in free space. However, the reason for the suppression with 
regard to the experiment of (Hood et al 2000) is actually the value of T a w IO7 at the 
antinode, thus D at w Df ree /10 Df ree . 




Figure 4. Diffusion coefficient D as a function of the axial position. Figure(a), 
parameters from (Hood et al 2000). Figure(b), parameters from (Pinkse et al 
2000). Same parameters as figure l2l From bottom to top: figure(a), D a t{cto, s), 
.D™.,, 116 el . Df ie i d + D a t(ao,s), exact (solid line), D a t {a, s ). Figure(b), D f ie id + 
fat(oo,s), Df ie i d + D at (a, s), exact (solid line), D„ n „, I16el . a is given by 11561 . 



5.2. Diffusion by eliminating the atomic variables 

There is also a possibility to include all quantities in the field variable, by setting 
a = in (|SJ| and inserting it in the equation for d. Similar arguments to the previous 
case lead to: 

H {Jield) = £ ^ c ( £C ) c t c ; K c (x)L c , a = c-a 
9 2 (x) 1 

' c a wl+^'l + s{x) (16c) 



g 2 {x) 1 



K c {x) = k + 7- 



w a + 7 2 1 + S ( X ) 

where a simplification relating a to s and p has already been made. The latter defines 
then the dipole force Fd = (3a) + (3* a. This case is trivial (we scale out the gradients): 

Df * eid = m2 afno c - {16d) 

We plot in figure^]!? = D fi e id+D a t(cto, s) for both experiments, in figurc^Jb) we also 
show D = D fi e id + D at (a,s). State (|15af) can also be used, the exact result will be 
approached differently. 

It should be clearly reminded that, in every case, Bloch equation-like and a 
damped displaced cavity mode are the background assumption. A picture is provided 
by imagining a swimmer swimming in a small pool, who consequently generates waves 
which get reflected by the finite size of the pool; that then changes the structure of 
the surface of the water to finally make the swim more difficult to the swimmer by 
pushing him according to the way he swims. In the end, the swimmer dissipates energy 
(r a 3> 7) in trying to move and gets tired. In the free space limit k — ► 00 , one has 
(f2 ,r ) — > (w ,7) and K c ~ k, and hence D at — » D free and D field — > 0. 
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5.3. Diffusion by considering the atom-mode coupled system 

Those previous expressions for the diffusion supposed the held or the atom as 
prescribed quantities. Other results are given where both objects are regarded as 
"active". This is done by solving JSJ) approximately, the procedure is analogous to 
(Hechenblaikner et al 1998). As the population difference a z is frozen, we consider 
two cases that could be grouped into one in future works. The result for the diffusion 
can be written D — D cav : 
with the "high saturation" case, 



\a \ 2 {l + s) 



or the "low saturation" case, 

_ |a| 2 |a| 2 gOQ nfiA 

^"^Hf (i + S )'hP(i + «) " 1 /j 

The difference is essentialy a factor 1/(1 + s), whose origin is related to the difference 
between < a^a > and < cr >< a >. An additive term j/\uj a \ 2 x |a| 2 /(l + s) has 
been dropped, it is a part of the free space expression and can be neglected within our 
approximations. If one sets s = 0, and explicits a by the first bounced state expression 
<|14a[l . then one recovers the expression derived by (Hechenblaikner et al 1998). The 
authors interpreted such term as cavity-induced, we explicit here its cooperative origin, 
indeed it is proportional to the imaginary part of v. For k — ► oo, $s(i/) — > and 
D cav cancels. Equation Qlfiefl is plotted in figure 01 Eventually, the diffusion due to 
spontaneous emission is proportional to the probability of being in the excited state, 
thus it is given by |9/|. 

6. Friction at first order in velocity 

The friction coefficient is derived from the same Liouvillians (|16of)(|16c|) as those for 
the diffusion, but we use another definition for the force. Since the friction coefficient 
is related to a velocity development, it is suitable to define here the force by minus 
the gradient of the Hamiltonians IlKjall (|l(jcf) . That force being given, we develop the 
expression to first order and relate its quantities to the spatial derivatives of steady 
state mean values. The method is a standard procedure, similar to (Gordon and 
Ashkin 1980, Hechenblaikner et al 1998). To perform the calculus, one suggestion is 
first to return from the variables (S, c) to (cr, a) in (jlfio| and (|16c|l . and then derive. 
The velocity v dependent force F v is given by the coefficient G (hereafter called 
friction) : 

F v = -(v-V 5 )V 5 G. (17a) 
6.1. Friction by eliminating the atomic variables 

The easiest case is the cavity- mode Liouvillian (|16c|) . which leads to a force 
proportional to a^a. By direct derivation, one obtains an important expression:* 

e c ( K c n c , , iu? c -k 2 c \ 

(176) 



G field - i\a\ ^ + R 2 [ n a +Jf a & + 2 ^ +R % Cc 
1 7 1 



£ c = 2g 2 " F(s, u), Cc = 2g ' F(s, v) , 

LU 2 + Y 1 + S UJ* + Y 1 + S 

Before the derivation, put (Q c — iK c )a fa E = constant 
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where s = 2g 2 \a\ 2 /\uj a \ 2 and (£ c (x),£ c (x) ,F(s,v)) are discussed below. For our 
regimes, it can generally be assumed Cc = 0. For (s = 0,F(s,u) = l), Gf ie id reduces to 
the friction obtained by (Hechenblaikner et al 1998) in the good cavity limit. 

6. 2. Friction by eliminating the field variables 

In the other limit i|16a|) . the friction obtained is easily interpreted as long as a free 
space term G/ ree is clearly recognized: 

Gat — Gf ree + £ a Gi + + (a(G 1( + £, a G2() 

U r 1 K 1 (17c) 

ga = 2g c F(s, v) , Ca = 2 g , 

+ 1 + S U}£ + K Z 1 + S 

where (Gf ree , G±, G2, Giq, G2Q) are functions of {su — 2\(3h\ 2 , T a , ft a , g) and given in 
|Appendix A.2| Similarly to Gfi e id, one can in general assume £ a = for G a t- For 
(£a = 0, Ca = 0), the free space term is recovered, but it is defined by the modified 
frequency f2 and decay rate T a . For (a — ► 0, — »■ 0, s = 0, F{s,v) — 1), G at /N 
(|A.5|I gives back the expression derived by (Hechenblaikner ef a/ 1998) in the bad 
cavity limit. 

In order totally to relate the friction coefficients (|17cllijl76(l to steady state 
quantities, one still has to specify the function F(s, v). 

6.3. Expressions for the function F(s, v) 

The dimensionless parameters £ a = — d g VL a , £ c = —d g Q c (and CqiCc) measure the 
sensitivity of the atom and cavity frequencies to the variations of the coupling. 
Some are oftentimes interpreted as trapping potentials U(x) acting on the atom, 
£ c = —dgU(x). In the low saturation limit, £ c reduces to 2gcj a /(oj a + j 2 ), and 
can be understood by analogy with a massive pointlike damped dipole oscillator 
(Hechenblaikner et al 1998). The inclusion of the saturation parameter, through the 
real function F(s,v)/(l + s), shows the departure from which the two level structure 
appears. The function, 

F(s,u) = l--±-X(gd g l3/0) , (18) 

-L — r~ S 

reduces to 1 at zero saturation, and is related to the derivative of (3. If the steady 
state is provided by the ping-pong schemes, then F(s, v) contains the variation of the 
referred states. However, in general, there is a possibility to provide several tractable 
and efficient expressions for F(s,u). An approximation is given which is part of a 
generalized ping-pong scheme to be dealt with elsewhere. One starts from a general 
bounced state expression, a — ao/(l — v p ), and a polarized state a — ao(l + v p ), with 
v p = vj{\ + s p ), and where s p (x) is the saturation parameter relative to a referred 
state "p". The derivatives (gd g /3)/fl become: 



2f, 



2v n 



gd g (3/(3 = 1 + -^F(s p , v) , gd g p/P = 1 + — ^F( Sp: v) , (19) 

bounce 1 — Vp pol. Y^rVp 

and hence F(s, v) (|18|) is related to the function F(s p , v) l|19|l relative to the referred 
state "p". The next step is to understand what kind of approximations are most 
efficient; different possibilities are given: 




Figure 5. Friction coefficient G as a function of the axial position. Figure(a), 
(Pinkse et al 2000). From bottom to top: Gyy ee (s), G a t{s), Gfieid, exact G (solid 
line). We used state 11561 and the v dependent version 12061 for the function 
F(s,v). Figure(b), (Hood et al 2000). From bottom to top: G fi e id + G a t(s), 
exact G (solid line), term £ a Gi + £^G2 + Ca(G 1( j + £ a G 2 £ ). We used state I15fel 
and F(s,u) is given by I20al . Notice, in figure(b), the positive contribution of 
the translational friction term. 



gdg/3/0 ~ 1 (20a) 
gd g (3/P ~ 1 + -21- 1 + _L_ (206) 

gd g (3/f3 ~ 1 + * 1 ~~~ (20c) 

1 — v a 1 + s 

gd g f3/(3 ~ 1 + J^-F(s,v) , (20d) 

where v a = vj(\ + s). The referred state p has been eliminated. The first case (|20a|l 
assumes a variation proportional to that of a pure standing wave. The atom could be 
affected by the reaction of the field but, when it moves, it sees a constant strength 
(/3 ~ g x constant) . Such approximation is generally valid for small values of v or at 
high saturation (check with l|19| l). as for (Hood et al 2000). In that case, H20af> gives 
F(s, v) — 1/(1 + s). The two approximations in l|206|l are obtained with the polarized 
state expression, it is assumed that the saturation parameter s p w s varies with v as 
s+is < 1, and F(s p , v) = 1/(1 + s) (pure standing wave). Such approximation is more 
valid at intermediate saturation v w s < 1. The arrow in equation l|206|l means that, 
at this step, v = 1 can be assumed. Indeed, the functions (£ a , £a> £ C) Cc) cancel for 
g = and hence the relevant points are around the maximum of the coupling. In such 
a case, ()18H simplifies to give F(s, v) = 1/(1 + s) 2 , and the difference of the exponent 
(l + s) with the function F(s, v) derived from I|20q|) shows again the limit between high 
and low saturation. The third case (|20c|) is aimed to work efficiently for large values 
of v, it is derived from the general bounced state expression. By the right arrow, 
we simplify the expression by taking the limit v p — > oo, and a supplementary ratio 
1/(1 + s) has been added from F(s p , v). The two arrowed limits just discussed have an 
interest (a part simplification) because they are smoothed versions of the expressions 
from which they derive, that might contain some undesired resonance. Finally, (|20dll 
represents the case close to the optical bistability way-of-thinking, the value of f3 is 
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Figure 6. Diffusion and friction coefficients as a function of the axial position. 
Figure(a), (7 = 0.0236, ft = 0.13, w a = 6, uj c = 1/6, go = 1, N = 2.16). From 
bottom to top at x = 0: D fi e id + D a t(s), exact D (solid line), D cav\16j\ , exact G 
(solid line), G fi e id- For F(s, u), we used the v dependent form of I20cl . Figurefb~). 
same parameters as figure 15Tb). From bottom to top at x = 0: Df ield + D a t(s), 
exact D (solid line), D e au (l6^ , Gfi^d, exact G (solid line). We used the v 
dependent form of 12061 . For both figures, we used state 115 at . 

returned, thus by H18|) . F(s, v) can be deduced (not written). Once these procedures 
made, the functions F(s, v) now depend on the saturation parameter of the steady 
state only, therefore one can simply forget the origin of the states that lead to <|19|) . 
As final simplifications, for the coefficients (G/ ree , G\, G%, Gi£, Gx) i|A.5|l . one can as 
usual operate Sh — s (as for the diffusion), thus all the friction coefficients are now 
completely determined by providing a. 

For the experiment of (Hood et al 2000), the translated terms of G a t (|17c|l have 
an effect and tend to render the friction positive. This is shown in figure [5] where we 
also plot the friction for the experiment of (Pinkse et al 2000). The friction Gfi e id 
\17b\ is an excellent approximation for that experiment. It shows that the quantum 
nature of the mode is again not needed, thus responding to the questioning of (Doherty 
et al 2000). A basic result for this paper is plotted in figure it shows a situation of 
high positive friction and a diffusion similar to (Pinkse et al 2000). The temperature 
is about twice as low. Also notice that, for figure GJa), it is Gfi e u that works despite 
the fact that 7 < k. This is a reminder that the bad/good cavity limits are only 
rough limits, in figureEJa) the opposite situation occurs. That is caused by the other 
physical parameters, like the coupling g, which are greater than k. Actually, this 
exact reversion can be explained by extending the notion of bad/good cavity limits: 
(T Q , K c ) are the relevant decay rates to determine which limit to use. In figure[lja), we 
have in fact T a > K c , thus the hierarchy is reversed (7 < k). Eventually, we present 
in figure a cavity resonant case, with a high photon number Nq~ 11, high friction, 
and a diffusion as low as (Hood et al 2000). In that case, the atom is trapped almost 
exactly half way between the node and the antinode, thus it experiences a non zero 
friction and diffusion (\7g ^ 0). Such parameters should be very efficient for trapping, 
with a very clear drop at the output transmission. 
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x [ W2 ] x [ X72 ] 

Figure 7. Diffusion and friction coefficients as a function of the axial position. 
Cavity resonant case. (7 = 0.18, k = 0.09, u> a = — 5, u) c = 0, go = lj ^0 = H)- From 
bottom to top at x = 0: Figure(a), D C av {16^ , exact D (solid line). Figure(b), 
G fi e id, exact G (solid line). For F(s,u), we used the v dependent form of 
120 cl . with state I14cl . The atom should be trapped at points xm such that 
<? 2 (xm) = with a temperature D/G ~ 2/35go — K. 

7. Remarks on optical bistability theory 

Our goal was first oriented towards cavity QED trapping mechanisms. It was only 
when v JJJl was labelled and attributed some importance, through the first ping-pong 
schemes (|13afll|136|) . that its interpretation became needed. Its link to the cooperative 
parameter C = g 2 /2n^ (Bonifacio and Lugiato 1978) is obvious and it is standard 
belief that v is a natural generalization of C. What is uncommon, however, is that v 
is never used in the optical bistablity theory (Lugiato and Narducci 1990, Abraham 
and Smith 1982), as far as we are concerned. We make essentially two remarks, 
showing that, maybe, one could rethink optical bistability with v and C in hand. 
There is a standard mean-field optical bistability state equation, which is commonly 
written in the form (Lugiato and Narducci 1990): 

x = y{l + i9 + 2C(1 - iS)/(l + 8 2 + M 2 )}- 1 , (21) 

where x, y are the scaled output and input field amplitudes respectively. That equation 
depends on the scaled detunings 9 = lu c /k, S = uj a /l an d the cooperative parameter 
C . Therefore, the space of parameters in which one works is the (C, 6, 8, y) space; 
however, as one parameter is usually a variable, like for example the input amplitude 
y, bistability conditions are drawn in the cubic space (C,8,9). We showed that the 
state equation above is actually the bounced state ©, scaled differently with the 
schematical correspondences y — > eta, x — > a b, and with a saturation parameter 
related to \x\ 2 . The referred state is here the state itself. Consequently, the space is 
now reduced to (v, ao) or (v) alone if ao is a variable. What was to be expected does 
happen; we show two typical inequalities that are conditional to optical bistability 
(Agrawal and Carmichael 1979, Drummond and Walls 1981): 

4[S9 + C-l] 3 >27C(5 2 + l){9 2 + 1), 2C > 59 - 1 . (22a) 

Those two conditions could be re-written, 

[M 2 + 25R(^)] 3 > 27M 4 , \v\ 2 > . (226) 
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That is a strong result. Those bistability conditions depend essentially on the value 
of v, thus we gain one dimension since there is a possibility to draw the conditions in 
the lA-complex plane. Notice incidentally that the conditions above relate the phase 
of v to its modulus. How far such simplifications go should be checked in order to 
determine the exact status of v and C. 

8. Conclusions 

If there was one conclusion to draw it would be the reminder that, even if quantum 
effects happen to be in the cavity QED setting, the confinement of the electromagnetic 
field within boundaries gets along with a strong classical behaviour. That can be 
interpreted as, due to the boundaries, even one photon on average would meet several 
times the atom, the latter sees a lot of photons, thus most likely a classical field. 
We showed that, even if one tunes the external laser quasi-resonantly with a dressed 
state, the coherent driving has a continuous component that almost washes away 
those quantum states. That is true for the steady state and for the diffusion and 
friction coefficients. We provided a self-consistent method that recognizes a new 
physical parameter and that allows one to extend the semiclassical limit out of both the 
bistability state equation domain and the low saturation regime. We derived classes 
of steady states that incremently took into account the quantum nature of the atom, 
while the field remained classical. Besides, we gave diffusion and friction coefficients 
that do not depend on any particular form of the steady state. We pointed out that 
the usual bad/good cavity limits can be reformulated in the strong coupling regime 
by the statement: The modified decay rates (T , K c ) become the relevant decay rates 
to determine which limit to use. 

The steady states are functions of the referred states, which in general are steady 
states but for other parameter regimes. The steady states belong to two distinct 
families, the bounced states, and the polarized states. They respectively refer to the 
saturation parameter and the polarization of the lastly referred state. The parameter 
v is a step that allows one to move from one state to another, it also represents the 
quantity that sends back the driving strength to the mode (or the atom) from E to 
vE. Thus, the method is called the Ping-Pong since the coupling g(ato~ + aa^) may 
represent the net of a ping-pong game between the atom and the mode, where (E, v) 
respectively represent the ball and the transporter. Although v does not depend 
explicitly on E, it however relies on a finite value of the driving strength £ / 0. 
Such relationship between v and E has as its cornerstone the lowest dressed state. 
The condition E — gets along with v = 1, the latter equality meaning that the 
lowest dressed state |— ) is tuned into resonance. In parallel, v seems to have a deeper 
meaning than the cooperativity parameter C. We showed that known inequalities that 
are conditional to optical bistability are actually functions of the complex parameter 
v only, whereas they are usually written in terms of (C, S, 9) . 

We want to notice that the pictures which appeared in this paper, whether they 
are the ping-pong or the swimmer, are descriptive of the structure of the system. They 
are believed to be complementary to the sisyphus effect which is more a dynamical 
picture of the atom climbing up and down potential hills. A structure effect has been 
cxplicited by showing that the reduction of the diffusion for the experiment of (Hood 
et al 2000) can be explained by a scaling argument on the dynamical variables due to 
the interactivity between the atom and the mode, and so such a reduction does not 
seem to be a mystery anymore. Similarly, the method is able to recover accurately 
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the experiment of (Pinkse et al 2000). We also derived a state, with a classical 
assumption for the field, which reproduces reasonably well a heterodyne transmission 
measured by (Hood et al 1998) whereas the bistability state equation is no longer 
valid. The understanding of the structure of the system should be deeply thought in 
future works because it lays the stress on settled (permanent) quantities, that are fixed 
as soon as dimensionless parameters are. We shall show elsewhere that there are an 
infinite number of steady states exhibiting fractal structures, which could hopefully 
account for instabilities in cavities. The ping-pong method is extendable to linear 
systems having a coupling term of the form discussed in this paper and thus a wide 
range of electrodynamical Liouvillians could be addressed. The classical assumption 
for the field is relevant for this paper only, the ping-pong method can be applied to 
quantum mode regimes. In our opinion, should future experiments in cavities be done 
to trap an atom by a quantum field, then one should test first a clear signature of the 
quantum origin of the mode by making for example a heterodyne and photodetection 
measurements. If those signals are clearly separated, the atom will possibly be trapped 
by a quantum field. None of those experiments clearly fulfilled such kind of conditions. 
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Appendix A. Results for the diffusion and friction coefficients 

Appendix A.l. Diffusion coefficient, atomic part 

Solution for the diffusion coefficient that is a generalization of (Gordon and Ashkin 
1980), for the standing wave case, with the Liouvillian given by: 

H = Q a (x)a j 'a - ga h (x)a f ~ ga* h (x)a , T a (x)L a , (A.l) 

where there is no need to specify the exact form of the frequency f2 a (x) and the 
decay rate T a {x). The force is given by l|166|l . The procedure is identical to (Gordon 
and Ashkin 1980), at the exception that the distinction between ah and ad leads to 
two suplementary terms here grouped in Dq. The result is for the diffusion l|16cf) . 
D = D at + 3?(£> ) : 

Do = k ng + 1* i + a s h ( 1 " (TT^ ( 1 + (n. 4 -h 2 + 2SA fi^Vj) ' (A ' 2) 

where A a = a\a* d 2 — |ctfe | 2 |«d | 2 . The free space like term is given by: 

= 5^,^(1 + (5$, - IK «4 + , (A.3) 

where Sh = 2\Ph\ 2 with j3h = gah/d a , ft a — ^ a ~ iT a . (|A.3(1 is sufficient with st = s 
and ad = {ol or ao). 
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Appendix A. 2. Friction coefficient, atomic part 

The friction is calculated by starting with (|16afl . or directly by IjA.ll) with — a . 
The force operator defined by minus the gradient of the Hamiltonian is developped to 
order one in velocity: 

Gat = -a Q al - a^a v - ^ a (x)(cr ] ' a) v , (A. 4) 

and is expressed in terms of the derivatives of the steady state values of (a, o~\ a 'a). 

{a^a) v = 7^—^ — ( - d g < a z > s +m\a*gd g < a > s /O a 



4r a 1 + s h 
a v = —2p h (a^a) v + id g < a > s /f2 a . 

Use I3 h = ga /& a , an d its derivative gd g (3 h = (3 h (l + g(£_ a + Ka)l& ah w here £ a , ( a are 
related to F(s, v) in (|17cll . The expression for the coefficient G a t l|17cf) finally reads: 

i 2| ao | 2 r»a r 2 , „,., , r 2 a 



a,M = al+T , fn _ p -(-«,, + 2 (i-. k )j 5Tif ) 

G2[sh) = ng+if (IW rli 2 s * + nfTff) (A - 5) 
r r w 1 2|an|2 g " a fi 4r ° 1 

G2dsh) = ' 

The analytical expressions for G a t,Gfi e id are explicited by providing a particular 
formula for F(s, v). The form of the functions tt a (x),T a (x) in (|16ajl only affect the 
functions £ a (%), (a(x). 
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